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QJ , Abstract. We develop the theory of interaction of classical plasma with Kaluza- Klein 

I (KK) gravitons in the ADD model of TeV-scale gravity. Plasma is described within the 

' kinetic approach as the system of charged particles and Maxwell field both confined on 

the brane. Interaction with multidimensional gravity living in the bulk with n compact 
extra dimensions is introduced within the linearized theory. The KK gravitons emission 
f~| ■ rates are computed taking into account plasma collective effects through the two-point 

correlation functions of the fluctuations of the plasma energy-momentum tensor. Apart 
^ i' from known mechanisms (such as bremsstrahlung and gravi-Primakoff effect) we find 

1^ . essentially collective channels such as the coalescence of plasma waves into gravitons 

which may be manifest in turbulent plasmas. Our results indicate that commonly used 
rates of the KK gravitons production in stars and supernovae may be underestimated. 



> 

in 

5^ ' 1. Introduction 

^ TeV-scale gravity proposal due to Arkani-Hamed, Dimopoulos and Dvali (ADD) [1] 

O ■ suggests that the standard model particles live in the four-dimensional subspace (the 
brane) of the D-dimensional bulk with n = D — 4 extra dimensions compactified on 
a torus which are inhabited only by gravity. The D-dimensional Planck mass 
is supposed to be TeV-scale, and the large extra dimensions (LED) to have sub- 
millimeter size. In this scenario there is an infinite tower of massive KK gravitons [2] 
whose existence may be detected using table-top and collider experiments, astrophysical 
and cosmological observations. It has been pointed out that one of the strongest 
bounds on the parameters come from the supernova SN1987A data [1]-[I0]. Processes 
contributing to energy loss due to graviton emission from SN1987A, the red giants and 
the Sun include photon-photon annihilation, electron-positron annihilation, nucleon- 
nucleon gravitational bremsstrahlung, gravi-Compton-Primakoff scattering [3]-|10j. 

Previous calculations of the KK graviton emission rates were performed using one- 
particle Feynman rules with subsequent averaging over Bose-Einstein and Maxwell- 
Boltzmann distributions to get the results reliable for finite temperature T [Ij. These 
calculations, however, did not take into account plasma collective effects such as 
Debye screening, interaction of charged particles with plasma waves etc, which may be 
important in astrophysical conditions. This is the goal of the present contribution. We 
generalize kinetic theory of interaction of plasmas with gravitational waves developed 
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earlier in 3 + 1 dimensions [TT] to the ADD model and present calculation of gravi- 
bremsstrahlung rate from non-relativistic thermal isotropic electron-ion plasma. 

2. Kinetic approach 

Consider coUisionless plasma consisting of charged particles of types a with parameters 

(=_ rr),_ HfiSPrihpH hv mirrnsmnip diRtrihiitinn fiinptinn II 9.\ II .31 



Cq, rua described by microscopic distribution function 

p) = ^ S{r - r,(t))5(p - p,(t)), (1) 



normalized as J J-'a{x,p)dp d^x = N^, satisfying the kinetic equation 

f^ + vf^ + e„(E+|vB])^^0. (2) 

They interact via the electromagnetic field E, B satisfying Maxwell equations 

divE = 47rp, rotE = divB = 0, rotB = 47rj + — , (3) 



(4) 



at ' ' 'at 

with the source terms 

p{x) = ^ea J^a{x,p)dp, j(x) = X^e^ / ^J^a{x,p)dp. 

To describe gravitational radiation we will need the 3-spatial components of the plasma 
energy momentum tensor 

where 

^'^ij = 5Z / ^-^«(^' P) ^P' Pa = Vp"^ + ^ = r}, (5) 

a J Pa 

pTij = {^i^j + BiBj - 5ij i ^ . (6) 

To solve the system of equations (IllHj) we use perturbation theory in terms of the 
electric charges. First we separate fluctuations from the average distribution using the 
approach of [HI [12] : 

J-(t,r,p) = /° + 5/° + 5/„, (7) 

where /° =< p) > - the equilibrium Maxwell distribution function: 

5/°- represents fluctuations due to chaotic particle motion ("zero" fluctuations), while 
5/Q,(t,r,p) - stands for fluctuations due to electromagnetic interaction of the particles. 
Using the fact that /° satisfies the free equation (for zero charges), we obtain: 

^ + v^ + ,„(E+lvBl)M±M±^«.o. (9) 

ot or op 
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We then further expand the fluctuation 6fa in power series in terms of charges: 
6 fa = + + and use the Fourier-transformation 

fix) = / f{k)e-''^^'d'k 

to get in the first two orders 

5/a(fc,P) = -(^;;eb)F(A:)^, F^E+[vB], (10) 

= ~ (27r)^(^-kv) / ~ P) + ^/<i(^i' P))^'^i- (11) 

The corresponding expansion of the charged particles energy-momentum tensor will 
read: 

oo „ 

mT^k = mTl + S^Tl + ^^/«(^' P) ^P- (12) 

First non-zero contribution to gravitational radiation comes from terms of the second 
order in the fields E, B, namely 6Tik = 5 mTf^ + 5 pTik- 

3. Gravitational radiation 

In ADD model one considers [2] the linearized D = A + n dimensional Einstein equations 

where ipMN = — VMNhp/2, Hmn = Qmn — Vmn, the indices M, N run over the 
brane n, u = 0, 1, 2, 3 and n directions on the torus, and the harmonic gauge d^il)^'^^ = 
is understood. Formula for the total (integrated over space and time) energy loss on 
gravitational radiation in ADD model was recently derived in 



where TsN{k) is the four-dimensional Fourier-transform, and 

j^SNLR ^ 1 y^SL^NR _^ ^SR^NLl^ _ _J_^SN^LR j^j^^^ 

is the polarization projection operator. Since the stress tensor has only brane indices 
jj, = 0,i and it satisfies the conservation equation k^T^'^{k) = 0, one can eliminate its 
components T^" in favor of T^^ . The resulting expression will contain the sum of various 
order terms from which the lowest non-zero contribution comes from the product of the 
second order terms. We also have to normalize the energy loss of the stationary plasma 
per unit volume and per unit time, which is achieved by introducing the two-point 
correlation function according to the relation 

< Tap{k)T:,p,{k + AA;) >=< T^^T^,^, >u 6\Ak). 
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The final expression for the gravitational emission rate to all graviton modes per unit 
plasma volume per unit time will read 



Nezd 



P = jTT^ E / < 6T,kST:,,, >u A"''"' , , (15) 



ji^iki k _ I \ ii \kk , \ik \ i k 

~ 2 



]^^^^^'^kk' ^^^k'^.'k-\^_^_^^k^^'k'^ A,, = - n,n,, n = -. (16) 

In our approach we did not use explicit decomposition of the full multidimensional metric 
perturbation Hmn in massless and massive modes, these correspond to N = and N 7^ 
terms in the sum respectively. For massless modes the three-dimensional vector n is the 
unit vector, but not for massive modes. Note that Hmn has bulk components due to 
the trace term in the wave equation. 

To calculate gravitational bremsstrahlung from the non-relativistic plasma we need 
the longitudinal field and the related fluctuation function. Keeping as the source in 
Maxwell equations terms of the zero and the first order 



divE = 47re j dw{5f + 5/^), 5f = 



e 1 ^ df 



Efc^, (17) 



muj — kv (9v 

one can express the first order electric field through zero fluctuations: 

47re2 [ , (kv) , 9/o(v) , , 

ELik) = 1 + --^ / dv ^ I 19 

muj\K\^ J (w — kv) c/v 

the last line being the longitudinal permittivity of the isotropic plasma. Using this, one 

then finds the two-point correlation functions defined above [13] Jil Correlation function 

for zero-order particle density < >° is calculable for the Maxwell distribution ([8]) in 

the finite form, while the field correlation functions can be then expressed through this 

quantity. Using the equations (fT7j) and (fTSll one find the following relations between 

Substituting this to (fluctuations of) the energy-momentum tensors and making use of 
the formula ( ITSl) . one obtains after some rearrangements: 

f,% U,,,^L1^ < ^s^,,,,)l^l^\ . (21) 

J [ m \eL{k2)\^ \^L{ki)r ^^(^ajr J 

According to the relations (!20|) , the first term here is the product of the field and particle 
density correlation functions, while the second is quadratic in the field correlators. The 
following frequency regions has to be distinguished: 

I Density fluctuations can be characterized by correlation functions whose Fourier-transforms exhibits 
homogeneity and isotropy: < Sn{ri,ti)6n{r2,t2) >=< 6n^ >r.t, r = ri — r2, t = <i — t2. The space- 
time Fourier transformation will give the spectral fluctuation density: < Sri^ >k— J J c?te^''"'+"^'' < 
6n^ >r,t, < 6n{k)6n{k') >= (27r)4 < Sn^ >k 5'^{k - k'). 
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i) Zero real part of the complex permittivity ReeL{u,k.) = in the denominators 
corresponds to propagation of the longitudinal Langmuir waves with the frequency 



V m 

and the long wavelength: |k| < where r^,^ = ^ is the Debye shielding radius. 

ii) For OJ > col, \k\vT^ ^ 1, then ~ 1 which corresponds to electromagnetic 
intaraction switched off. 

iii) For u < ul one has ^^(fc) — 1 + (|k|r£i)^^ which corresponds to the Debye 
shielding. 

Correspondingly, we have to distinguish radiation losses due to gravitational 
bremsstrahlung (denominators are non-zero, both terms being taken into account), and 
the coalescence of two Langmuir plasmons into graviton ( the second term with both 
denominators at Langmuir frequencies). 

3.1. Gravitational bremsstrahlung 

This corresponds to Reei 0. The main contribution to the integral comes from the 
region u < \\<.i\vt^, so for > ul, |ki|r£) > 1, then — 1. In opposite limit /cq < wl, 
Debye shielding becomes manifest, then 5l(^i) — 1 + (|ki|r£))^^. Therefore we obtain 
in the two cases: 

P^-^(l + v^)eXV^{/i(n)lnA + -J2H}, > co,, (22) 

1 oo . .n+l 

P^-|(l + v^)eXV^/i(n)lnA, k' < u,, (23) 



where In A is the Coulomb logarithm, A = r£,^/r, r^^ = ^ - is the electron Debye 
radius, f is determined from the condition e^/f = T and M^, = ^ is the D- 
dimensional Planck mass. The coefficients Ii{n), l2{n) depend on the number of extra 
dimensions n, 2 < n < 7 as follows: {h{n) = 1.1, 1.07 0.97 0.83, 0.67, hin) = 
4.8, 6.35, 7.3, 7.5, 7.1. If n = {Ii{n) = 1, hin) = 1) these results coincide with 
those - in 3 + 1 dimensions [11]. Without collective effects, the result in 3 + 1 dimensions 
corresponds to the quadruple formula calculation with consequent thermal averaging 

3.2. Coalescence of two Langmuir plasmons into graviton 

For coalescence of two Langmuir plasmons into the graviton the field correlation function 
can be presented as 

<E^ >k=A'K''T{5{uj-UL) + S{u + UL)) (24) 

The graviton frequency will be twice the Langmuir frequency and the rate of emission 
will be 

^° = 2e^^eV^^^/., (25) 
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with J„ = 1.4, 1.9, 2.2, 2.3, 2.2, forn = 2, . . . , 6. 



Compare with the resuh of the ref. [4] for the bremsstrahlung losses which was 
obtained using the one-particle approach: 

^ r(| + n)n,n,ZyiGB{n) T^+^ 

Hi) ^ 



3 



^2 



2 



where n^. and are the electrons and ions densities, Zj is the ions charge, a = e 
and the numerical integral Icsin) is presented for n = 2, 3 : IcBin) = 0.7, 0.3. Being 
applied to red giants with the electron temperature T ~ 8.6 keV and the electron density 
rig = 3.0 ■ lO^^cm"^, assuming the conservative upper limit on the energy-loss rates of 
red giants e ~ lOOer^f g~^ sec~^ , the authors of [1] obtained the upper limit for 
Planck mass to be 10^ for n = 2 and 1.7 ■ 10~^ for n = 3. Our result fl23|) for the 
bremsstrahlung losses is about three time greater than ( l26l) . Moreover, an additional 
enhancement comes from the in coalescence process (125|) . 
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